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In the spacetime of a d-dimensional static tense brane black hole we elaborate the mechanism 
by which massive scalar fields decay. The metric of a six-dimensional black hole pierced by a 
topological defect is especially interesting. It corresponds to a black hole residing on a tensional 
3-brane embedded in a six-dimensional spacetime, and this solution has gained importance due to 
QQ ] the planned accelerator experiments. It happened that the intermediate asymptotic behaviour of the 

, fields in question was determined by an oscillatory inverse power-law. We confirm our investigations 

' by numerical calculations for five- and six-dimensional cases. It turned out that the greater the 

I brane tension is, the faster massive scalar fields decay in the considered spacetimes. 

o 

Q 

■ I. INTRODUCTION 

In recent years the subject of the late-time behaviour of various fields in the spacetime of a collapsing body has 
acquired a great deal of attention due to the fact that regardless of details of the collapse or the structure and 
O ^ properties of the collapsing body the resultant black hole may be described by a few parameters such as mass, charge 
D I and angular momentum, hlack holes have no hair. 
j-^' In |1] it has been pointed out that the late-time behavior is dominated by the factor ^-(2^+3)^ foj. each multipole 
moment I, while in that the decay-rate along null infinity and along the future event horizon is governed by the power 
^— I i: laws and v^'^^^'^\ where u and v were the outgoing Eddington-Finkelstein (ED) and ingoing ED coordinates. 

I> ' The scalar perturbations in the Reissner-Nordtsrom j^^) spacetime were treated in [31 . It happened that a charged 
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PACS numbers: 



, hair decays slower than a neutral one Q-p- Burko [3] reported the late-time tails in gravitational collapse of a self- 
^ ■ interacting (SI) fields in the background of Schwarzschild solution. On the other hand, RN solution at intermediate 
, late-time was considered in Q . The nearly extreme case of the RN spacetime was elaborated analytically in article 
Q.It was concluded that the inverse power law behavior of the dominant asymptotic tail is of the form t^^/^ sin(mi). 
■ It turned out TO] that the oscillatory tail of scalar field in Schwarzschild spacetime has the decay rate of t~^'^ at 
asymptotically late time. The power-law tails in the evolution of a charged massless scalar field around a fixed 
background of dilaton black hole was studied in jTl| , while the case of a self-interacting scalar field was elaborated 
in |12l | . The analytical proof of the above mentioned behaviour of massive scalar hair in the background of a dilaton 
black hole in the theory with arbitrary coupling constant between U{1) gauge field and dilaton field was given in 
[isj . It was revealed [3| that the asymptotic behaviour of massive Dirac fields in the background of Schwarzschild 
; ^ ' black hole was independent of the multipole number of the wave mode and of the mass of the Dirac field. Their 
Ci [ decay was slower comparing to the decay of a massive scalar field. The behaviour of Dirac fields in the spacetime 
of RN black hole was treated in [l5'| while the case of a stationary axisymmetric black hole background was studied 
numerically in . In [17| both the intermediate late-time tail and the asymptotic behaviour of the charged massive 
Dirac fields in the background of a Kerr-Newmann black hole was under consideration and it was demonstrated that 
the intermediate late-time behaviour of the fields in question was dominated by an inverse power-law decaying tail 
without any oscillation. Massive vector field obeying the Proca equation of motion in the background of Schwarzschild 
black hole was studied in jl^ where it was revealed that at intermediate late times, three functions characterizing 
the field have different decay law depending on the multipole number I. On the contrary, the late-time behaviour is 
independent on i.e., the late-time decay law is proportional to t~^/^ sin(?7it). In [lH the analytical studies concerning 
the intermediate and late-time decay pattern of massive Dirac hair on the dilaton black hole were conducted. Dilaton 
black hole constitutes a static spherically symmetric solution of the theory being the low-energy limit of the string 
theory with arbitrary coupling constant a. In [23| it was assumed that the massive Dirac hair propagated in a static 
spherically symmetric spacetime with asymptotically fiat metric. The intermediate late time behaviour was found 
numerically and the final decay pattern at very late times was calculated analytically. The metric in question, was 
characterized by the ADM mass and some other parameters of the background fields. The intermediate and late-time 
decay of massive scalar hair on static brane black holes was elaborated in , while the decay of massive Dirac hair 
in the spacetime of black hole in question was studied in [l^] . 
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For n-dimensional static black holes, the no-hair theorem is quite well established [235 . One should also mention 
about some attempts to establish uniqueness theorem for another kind of black objects, black rings. The uniqueness 
theorem for five-dimensional stationary the so-called Pomeransky-Sen'kov black ring was found in [2J] while black ring 
CT-model uniqueness theorem was given in [25j. On the other hand, the mechanism of decaying black hole hair in higher 
dimensional static black hole case concerning the evolution of massless scalar field in the n-dimensional Schwarzshild 
spacetime was determined in [2^ . The late-time tails of massive scalar fields in the spacetime of n-dimensional static 
charged black hole was elaborated in [27'| and it was found that the intermediate asymptotic behaviour of massive 
scalar field had the form ^-('+"72-1/2)^ This pattern of decay was checked numerically for n — 5,6. Quasi-normal 
modes for massless Dirac field in the background of n-dimensional Schwarzschild black hole were studied in |28| . 

Due to the fact that forthcoming Large Hadron Collider (LHC) may open the way for producing microscopic black 
holes which will quickly decay through the Hawking radiation, interest in tensional brane black holes rapidly grows. 
In [2^ the exact metric describing a black hole localized on a codimensional-2 brane was constructed. It turned out 
that the line element (Pild-2 on the unit (d — 2)-dimensional sphere is replaced by the metric on a unit sphere but 
with the wedge removed from the polar coordinates, i.e., one gets some kind of a topological defect piercing the black 
hole in question. The finite brane tension modifies the standard results achieved in the case of a brane black hole 
with negligible tension. 

This discovery bears directly on the question of the underlying physics. Namely, for the first time it was observed in 
the semi-classical description of the black hole decay process [30] . In [sij the metric of codimensional brane black hole 
was generalized to the rotating case, while article j32j was devoted to massless fermion excitations on a tensional three- 
brane embedded in six-dimensional spacetime. Various aspects of the tense brane black hole physics were studied. 
Among all, quasi-normal modes, bulk scalar emission, grey-body factors and the behaviour of scalar perturbations in 
the background of a black hole localized on a tensional three-brane in a world with two large extra dimensions were 
elaborated in [33|-(36|. Having all in mind it will not be amiss to pay more attention to the decay process of massive 
scalar hair in the background of a tense brane black hole. 

In this paper our investigations will be devoted to studies of d-dimensional tensional brane black hole. In the case 
of d = 5, 6 we provide numerical analysis of the intermediate late-time behaviour of massive scalar fields in the black 
hole backgrounds in question. 

The remainder of the paper is as follows. In Sect. [H] we shall discuss the massive scalar field behaviour in the 
spacetime of d-dimensional static spherically symmetric tense brane black hole. In Sect. IIIII we treat numerically the 
problem of five and six-dimensional black hole in question. We conclude our investigations in Sect. IIVI 



II. MASSIVE SCALAR FIELDS IN A TENSE BRANE BLACK HOLE BACKGROUND 

Recent progress in string theory and brane world models led us to an intense interest in higher-dimensional black 
holes. In model with large extra dimensions black hole size is much smaller than the effective size of the extra 
dimensions. Simply additional dimensions are considered as having infinite extent. Models with large extra dimensions 
acquire much more attention also because of an interesting possibility of lowering the fundamental scale of gravity 
down to order of TeV. It can be argued that future experiments, which involve high-energy particles collisions at 
future colliders, hold great promise for illuminating the nature of mini black holes. Especially, forthcoming LHC 
brings the possibility of creating microscopic black holes which will disappear quickly after creation with the emission 
of the Hawking radiation. Most examinations of extra-dimensions black holes and their physics were devoted to the 
zero brane tension case. But in principle finite brane tension may modify the physics of a tensional brane black 
hole. The nonzero tension on the brane can curve the brane as well as the bulk. Therefore more attention should be 
directed to the studies of the afore mentioned black objects. A tense brane black hole is locally a higher-dimensional 
Schwarzschild solution [2^ threaded by a tensional brane. This caused that a deficit angle appeared in the {d — 2)- 
dimensional unit sphere line element. In what follows, our considerations will be devoted to studies of such kind of 
black objects pierced by a codimension-2 brane. 

In our considerations we take into account the following equation of motion for massive scalar field: 

V*V,V^ - m^^ = 0, (1) 

where m is the mass of the field. One should study its behaviour in the background of static spherically symmetric 
d-dimensional black hole. The line element of such a black hole is subject to the relation 

ds' = -nr)dt' + ^+r'dnl^„ (2) 
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where / (r) = 1 — yro/T j , '"o is the radius of the black hole event horizon, while <ii^^_2 is a line element on 5*' 
sphere provided by the relation 



d-2 



d-3 i-1 
i=2 j = l 



(3) 



In our considerations we take into account S'^~'^ sphere threaded by a codimension-2 brane, so the range of one of the 
angles, let us say 0^, will he < (fn < 2itB. 

Consequently, in the case of d-dimensions we resolve the scalar field function as follows: 



(4) 



where YJ" are scalar spherical harmonics on the unit {d — 2)-sphere. For the case when (d — 2) sphere is pierced by 
a topological defect, one can introduce hyper- spherical coordinates in the form: < r < oo, < < n, — tt < 4>i < 
TT, . . ., —ttB < <j)d-3 < T^B. Having in mind results of (37j . one may introduce k = XliLo where rii are the integers 
separation constants and now the multipolc number I has the form 



I = k + n, 



(5) 



where n = ^ expresses the existence of a topological defect. Further on, it is convenient to make the change of 
variables. Namely, let us define the tortoise coordinates y 



dy = 



dr 



1 



d-3 



(6) 



Using Eq.([5]) and relation ^ one has the following equations for each multipole moment: 

ip,tt - ip,vv + Vip ^0, 
where the form of the strictly radial potential V{r) yields 



V = fir) 



d-2 



1 



f 2 dv 



r---' fir) + 



k + d-3- 



The square of the Laplace operator on S"^ ^ threaded by a topological defect is given by 

= { k+\fi\ 1 { k + d-3- 



(7) 



(8) 



(9) 



One can introduce an auxiliary variable ^ which satisfies ^ — fir)ip. Next, one assumes that the observer and the 



initial data are in the region where the following condition is satisfied r <^ 
intermediate asymptotic behaviour of massive scalar field, i.e., r <^ t <^ 



(romp 



and moreover we consider the 



7— — Yj . If we consider further, the explicit 
dependence on time in function ip of the form Q-i-omegat ^ where uj is frequency, then we achieve 



2 2 viv + ^) 
Ll> — m 



(10) 



where v 



Solution of Eq. PU)) can be found by virtue of the method used in [38] . It happened 



that the asymptotic intermediate behaviour of scalar field at fixed radius implies 



(11) 



while the intermediate behaviour of massive scalar fields at the future black hole horizon is dominated by an 
oscillatory power law tails of the form as follows: 



iP ^ sin(mi;). 



(12) 
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where v is is advanced null coordinate. 

In what follows we shall treat the problem of five and six-dimensional static spherically symmetric black hole threaded 
by a topological defect. The six-dimensional case is of a special interest. In spherical coordinates its line element 
implies 



dr' 



r^{de^ 



+ sin^ (j){dx^ -f sin^ x^V)])- 



(13) 



Now, the parameter B = \ 



27tM, 



measures the deficit angle about an axis which is parallel to the three-brane 



with finite tension A. Af* is the fundamental mass scale of six-dimensional gravity. The black hole horizon is at the 
distance of rg = Vg/B^^^, where Vs is the six-dimensional Schwarzschild radius given in terms of the ADM mass by 
the relation 



1 

4^ 



J_ (MbhY 



(14) 



One can remark that the main effect of the brane tension is to change the relation between the black hole mass and 
the event horizon radius. 

In the next section we shall check our predictions numerically confining our attention to the case oi d — b and 
d = 6. 



III. NUMERICAL RESULTS 



Due to the difficulties in solving differential Eqs. of motion in terms of adequate special functions for arbitrary 
spacetime dimension, we treat them numerically. By the method presented in [1, [H we shall find the solution of 
the equation of motion transformed into retarded and future null coordinates (u,v). Consequently, we achieve the 
following: 

4^-,™ + = 0. (15) 

Equation (I15p will be solved on an uniformly spaced grid using the explicit difference scheme. As was pointed out 
previously, the late time evolution of a massive scalar field is independent of the form of the initial data. In order to 
perform our calculations we start with a Gaussian pulse of the form 

^(u = 0,f)=Aexp(^-^^^— (16) 

Because of the linearity of the relation ([7]) one has freedom in choosing the value of the amplitude A. For our purpose 
we fix it as A = 1. The rest of the initial field profile parameters will be taken as vq = 50 and a — 2. 

We begin our numerical studies with the evolution of scalar field -0 on a future timelike infinity i+. In the consider- 
ations we approximate this case by the field at fixed radius y = 50. The numerical results for fc = 0, n = 1, B = 0.9 
and mass of the field m = 0.01 for different spacetime dimensions are depicted in Fig. [TJ One should remark that 
the initial evolution of the scalar field in question is determined by prompt contribution and quasi-normal ringing. 
Then, with the passage of time a definite oscillatory power-law fall-off is manifest. We get the following power-law 
exponents —3.14 and —3.64 for d = 5, 6, respectively. Having in mind Eq. ljlip one can remark that we get good 
agreement between numerical calculations and analytically predicted values: —3.11 and —3.61, accordingly, that is 
the discrepancies are only 1% and 0.8%. For all considered cases, the discrepancies between analytical and numerical 
values range from 0.2% up to 1.8%. The period of oscillation is T = tt/to ~ 314.5 ± 0.5. 

Next, in Fig. 2 we elaborate the evolution of the massive scalar field on the black hole future event horizon 
approximated by ip^u = 10'', t). We studied the evolution as a function of v for different spacetime dimensions, i.e., 
d = 5, 6. The calculations parameters were as follows: k = Q, n = 1, B — 0.9, m = 0.01. The power-law exponents 
and the period of oscillation are nearly the same as in Fig. [1] Then, we studied the behaviour of -0 due to the change 
of brane tension B. In our calculations we take B = 0.9, 0.7, 0.5 (for fc = 0, n = 1), while d — 5, 6, respectively. The 
obtained results are presented in Figs. [3]and|4l The received power-law exponents are —3.14, — 3.47, — 4.04 for the 
five-dimensional case and —3.64, — 3.97, — 4.54 for the six-dimensional spacetime. We observe that the bigger B is 
the smaller power-law exponents we get, i.e., the greater brane tension is exerted the faster decay of massive scalar 
hair on the considered black hole spacetimes is observed. 
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We also studied the power-law decay behaviour for different n = 0, 1, 2 (Fig. 5) and established B = 0.9 and A; = 
in five-dimensional black hole spacetime, where we obtained the power-law exponents equal to —2.03, — 3.14, — 4.26. 
It turned out, that the same calculation parameters reveal in six-dimensional case the power-law exponents of the 
form —2.53, — 3.64, — 4.78 (see Fig. 6). Figs. [7] and [5] depict the dependence of the decay pattern on the multipole 
number k. The case of different masses of the scalar field in the black hole background was treated in Figs. [9]and[T0l 
Namely, we studied the late-time decay rate on future timelike infinity in five-dimensional black hole spacetime for 
different masses of the field in question (see Fig. [5]). The decay rate was a power-law fall-off with the slope —3.14. 
The six-dimensional case was presented in Fig. 1101 For various masses of scalar field one gets the power-law fall-off 
with the slope of the curve equal to —3.64. 



IV. CONCLUSIONS 



Motivated by the forthcoming possibility of producing in LHC experiments microscopic black holes and the necessity 
of studying their properties we have elaborated the intermediate behaviour of massive scalar fields in the background 
of a static d-dimensional tensional brane black hole. Due to the tremendous difficulties in solving equations of motion 
for the massive scalar field in general d-dimensional black hole spacetime, we studied numerically the case of a five 
and six-dimensional spacetime. The case of a six-dimensional black hole pierced by a topological defect is of a special 
interest due to the interpretation of this line element as a black hole residing on a tensional 3-brane embedded in a 
six-dimensional spacetime. We checked analytical calculations by numerical ones and obtained good agreement with 
analytically predicted values. Numerical studies of the intermediate late-time behaviour of massive scalar field on 
future timelike infinity and on the future event horizon of tense brane black hole revealed the inverse power-law decay 
of the field in question. The inverse power-law decay pattern was checked numerically for d — 5, 6. It turned out that 
the intermediate late time behaviour depended on the mass of the scalar field, as well as the multipole number of the 
wave mode and parameters characterizing topological defect. The same situation takes place in the five-dimensional 
case. As we have mentioned, for a six-dimensional codimension-2 brane black hole its event horizon radius is strictly 
connected with three-brane tension. Thus, our numerical calculations revealed that the greater brane tension is the 
faster massive scalar field decays in the spacetime in question. 

To finish with let us give some remarks concerning the four-dimensional case. In four-dimensions the considered 
line element reduce to the case found in 39| . It describes a thin cosmic string passing through a Schwarzschild black 
hole. It was revealed that it constituted the limit of a much more realistic situation when a black hole was pierced by a 
Nielsen-Olesen vortex (40| . In the extremal limit of the black hole under consideration a completely new phenomenon 
occurred. The flux of the vortex is expelled from the black hole (for some range of black hole parameters), rather 
like the flux is expelled from a superconductor (the so-called Meissner effect) . It happened that in the case of the 
extremal black holes in dilaton gravity one has always expulsion of the Higgs fields from their interiors [4l| . 

Using the method presented in [sj , it could be easily found that the eigenvalues of the Laplace operator on a 
unit sphere pierced by a cosmic string, i.e., = {k+ \ h |) {k+ \ h \ +1), where fc = 0, 1, 2 . . . It turned out that the 
crucial role is played by the factor B connected with a mass per unit length of a string, because of the fact that 
h = n/B. Having in mind [l^, [l^l it can be shown that the intermediate asymptotic behaviour of a scalar field 
depends on the field's parameter mass, cosmic string parameter as well as the multipole number k. The late-time 
behaviour reveals that the power law decay rate which is independent of the parameters in question has the form of 
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FIG. 1: Evolution of the scalar field | | on the future timelike infinity i+ as a function of t for different spacetime dimensions 
d = 5, 6. Calculation parameters: m = 0.01, fe = 0, n = 1, B = 0.9. The power-law exponents are: —3.14 and —3.64 for 
d = 5, 6, respectively. The period of oscillation is T = — ~ 314.5 ± 0.5 for both curves. 
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FIG. 2: Evolution of the scalar field | V I on the black hole future event horizon H+ (approximated by = 10*, v)) as 
a function of v for different spacetime dimensions d = 5, 6, upper and lower curve, respectively. Calculation parameters: 
m = 0.01, fc = 0, n = 1, B — 0.9. The power-law exponents arc: —3.15 and —3.67, for the worst bottom curve. The period of 
oscillation is T = — ~ 314.5 ± 0.5 to within 1% for both curves. 
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FIG. 3: Evolution of the scalar field on the future timelike infinity i+ as a function of t for different B = 0.9, 0.7, 0.5 (curves 
from the top, respectively) and d = 5. Other calculation parameters: m = 0.01, fe = 0, n = 1. The power-law exponents are: 
-3.14, - 3.47 and -4.04. The period of oscillation is T = ^ ~ 314.5 ± 0.5 to within 0.2% for all curves. 
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FIG. 4: Evolution of the scalar field | V I on the future timelike infinity i+ as a function of t for different B = 0.9, 0.7, 0.5 
(curves from the top, respectively) in d = 6 spacetime. Other calculation parameters: m = 0.01, k = 0, n = 1. The power-law 
exponents are: —3.64, — 3.97 and —4.54. The period of oscillation is T = — ~ 314.5 ± 0.5 to within 0.4% for all curves. 



12 




FIG. 5: Evolution of the scalar field \tp\ on the future timelike infinity i+ as a function of t for different n = 0, 1, 2 (curves 
from the top, respectively) in d = 5 spacetime. Other calculation parameters are: m = 0.01, k = 0, n = 1. The power-law 
exponents are: —2.03, — 3.14 and —4.26. The period of oscillation is T = — ~ 314.5 ± 0.5 to within 0.3% for all curves. 
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FIG. 6: Evolution of the scalar field | V' | on the future timelike infinity i+ as a function of t for different n = 0, 1, 2 (curves from 
the top, respectively) in d = 6 spacetime. Other calculation parameters in Fig. 5. The power-law exponents are: —2.53, — 3.64 
and -4.78. The period of oscillation is T = -2^ ~ 314.5 ± 0.5 to within 0.6% for all curves. 
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FIG. 7: Evolution of the scalar field | V I on the future timelike infinity i+ as a function of t for different fe = 0, 1, 2 (curves 
from the top, respectively) in d = 5 spacetime. Other calculation parameters are: m = 0.01, n = 1, B = 0.9. The power-law 
exponents are: —3.14, — 4.15 and —5.20. The period of oscillation is T = — ~ 314.5 ± 0.5 to within 1.7% for all curves. 
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FIG. 8: Evolution of the scalar field | '0 | on the future timelike infinity i+ as a function of t for different fc = 0, 1,2 (curves from 
the top, respectively) in d = 6 spacetime. Other calculation parameters in Fig. 7. The power-law exponents are: —3.64, —4.67 
and —5.60. The period of oscillation is T = — ~ 314.5 ± 0.5 to within 3.5% for all curves. 
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FIG. 9: Late-time decay rate of the field | ip \max at i+ for the five-dimensional case and different masses of the scalar field m 
(values written above the lines). Only maxima of the oscillations are showed. The dashed line has the slope equal to —3.14. 
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FIG. 10: Late-time decay rate of the field | ip \max at i+ for the six-dimensional case and different masses of the scalar field m 
(values written above the lines). Only maxima of the oscillations are showed. The dashed line has the slope equal to —3.64. 



